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Abstract 

A second order, 0{a 2 s ), calculation in perturbative quantum chromo dynamics 
of the two particle inclusive cross section is presented for the reaction p + 
p —* 7 + c + X for large values of the transverse momentum of the prompt 
photon and charm quark. The combination of analytic and Monte Carlo 
integration methods used here to perform phase-space integrations facilitates 
imposition of photon isolation restrictions and other selections of relevance in 
experiments. Differential distributions are provided for various observables. 
Positive correlations in rapidity are predicted. 
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I. INTRODUCTION 



More precise examination of the expectations of perturbative quantum chromodynam- 
ics(QCD), including dynamical correlations inherent in the hard-scattering matrix elements, 
is made possible by the observation of inclusive production of two particles or jets each car- 
rying a large value of transverse momentum. Because they couple in a point-like fashion to 
quarks, the observation of photons with large values of transverse momentum in a high en- 
ergy hadron collision has long been regarded as an incisive probe of short distance dynamics. 
In the case of inclusive production of heavy quarks, the large mass of the quark and/or the 
fact that the quark carries large transverse momentum justifies use of a perturbative short- 
distance approach. In this paper, we continue our examination of the associated production 
of a prompt photon along with a heavy quark Data are beginning to become available 
on the associated production of a photon 7 carrying large transverse momentum along with 
a charm quark c whose transverse momentum balances a substantial portion of that of the 
photon [Q. An intriguing possibility is that the data may be used to measure the charm 
quark density in the nucleon. 

In this paper we report results of a full next-to-leading order perturbative QCD calcu- 
lation of p + p^'-f + c + X at high energy. For values of the transverse momentum p 1 ^ of 
the charm quark much larger than the mass m c of the quark, only one direct hard scattering 
subprocess contributes in leading order: the quark-gluon Compton subprocess gc — > 7c. The 
initial charm quark and the initial gluon are constituents of the initial hadrons. In addition, 
there is a leading order fragmentation process in which the photon is produced from quark 
or gluon fragmentation, e.g., gg —>■ cc followed by c — » 7X, or qc — > qc followed by q —>■ 7. 
At next-to- leading order in QCD, several subprocesses contribute to the 7 + c + X final 
state: gc — > gey, gg — > ccj, qq — > ccj, qc — > qcy, qc — > qcy, cc — > ccj, and cc — > ccj. A 
full next-to-leading order calculation requires the computation of the hard-scattering matrix 
elements for these two-to-three particle production processes as well as the one- loop 0{a 2 s ) 
corrections to the lowest order subprocess gc — > 7c. 
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We are interested ultimately in the fully differential two-particle inclusive cross section, 
E 1 E c da /cPpjCppc, where (E,p) represents the four- vector momentum of the 7 or c quark. 
For each contributing subprocess, this calculation requires integration over the momentum of 
the unobserved final parton in the two-to-three particle subprocesses (g, c, q, or q). Collinear 
singularities must be handled analytically by dimensional regularization and absorbed into 
parton momentum densities or fragmentation functions. In the theoretical analysis reported 
here, a combination of analytic and Monte Carlo integration methods is used to perform 
phase-space integrations over unobserved final-state partons and the momenta of the initial 
partons. This approach facilitates imposition of photon isolation restrictions and other se- 
lections of relevance in experiments. We work in the massless approximation, m c = 0. To 
warrant use of perturbation theory and the massless approximation, we limit our consider- 
ations to values of transverse momenta of the photon and charm quark p^ c > 10 GeV. 

In the lowest order direct subprocess, gc — > 7c, the prompt photon emerges in isolation 
from the only other particle in the hard scattering, the charm quark. Long-distance quark- 
to-photon and gluon-to-photon fragmentation processes have been emphasized theoretically 
H and parametrized phenomenologically in leading order || , and evolved in next-to- leading 
order These terms may account for more than half of the calculated inclusive single 

photon cross section at modest values of transverse momentum at the Fermilab Tevatron col- 
lider. Photons originating through fragmentation are likely to emerge in the neighborhood 
of associated hadrons. An experimental isolation restriction is needed before a clean iden- 
tification can be made of the photon and a measurement made of its momentum. Isolation 
reduces the size of the observed fragmentation contribution. Photon isolation complicates 
the theoretical interpretation of results, however, since it threatens to upset the cancellation 
of infra-red divergences in perturbation theory [7). In this paper, we include the fragmenta- 
tion contributions, and we impose isolation requirements through our Monte Carlo method. 

A combination of analytic and Monte Carlo methods similar to that we employ in this 
paper has been used to carry out next-to-leading order QCD calculations of other processes 
including inclusive prompt photon |§ and photon pair || production in hadron collisions, 
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single |T0] and pair production of heavy gauge bosons [11] , and heavy flavor pair production 



HI- 

In Section II, we describe the combination of analytic and Monte Carlo methods we use to 
carry out the next-to-leading order calculation. The next-to-leading order calculation itself is 
presented in Section III. Differential cross sections and other numerical results are discussed 
in Section IV. Summary remarks are collected in Section V. An appendix is included in 
which we derive analytic expressions for some of the parton level cross sections we use. 



II. MONTE CARLO METHOD 

The combination of analytic and Monte Carlo techniques used here to perform the phase 
space integrals is documented and described in detail elsewhere , so our discussion will 

be fairly brief, highlighting features important to our calculation. For the two-to-three parti- 
cle hard-scattering subprocesses, the technique consists in identifying those regions of phase 
space where soft and collinear singularities occur and integrating over them analytically in 
4 — 2e dimensions. In this way the singularities are exposed as poles in e. These regions are 
isolated from the rest of the three-body phase space by the imposition of arbitrary bound- 
aries through the introduction of cutoff parameters, 5 S and 5 C . The soft gluon region of 
phase space is defined to be the region in which the gluon energy, in a specified reference 
frame, usually the subprocess rest frame, falls below a certain threshold, 5 s V§/2, where 5 S is 
the cutoff parameter, and s is the center-of-mass energy in the initial parton-parton system. 
Labelling the momenta for the general three-body process by p\ + p 2 — > £>3 + £>4 + Ps, we 
define the general invariants by s^- = (pi + Pj) 2 and £y = (pi —pj) 2 . The collinear region is 
defined as the region in which the value of an invariant falls below the value 5 c s. 

The full set of three-body final-state subprocesses is: 

g + c— >g + c + 7 (2.1a) 

g + g -»• c + c + 7 (2.1b) 

g + g-*c + c + 7 (2.1c) 
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g + c^g + c + 7 (2. Id) 

g + c^g + c + 7 (2.1e) 

c + c^c + c + 7 (2. If) 

c + c^c + c + 7 (2.1g) 

The matrix elements integrated over the mutually exclusive soft and collinear regions 
of phase space are not the full two-to-three body matrix elements but, instead, specific 
approximate versions. In the soft-gluon case the approximate version is obtained by set- 
ting the gluon energy to zero everywhere it occurs in the matrix elements, except in the 
denominators. For the collinear singularities, each invariant that vanishes is in turn set to 
zero everywhere except in the denominator. This form is the leading pole approximation. 
The phase space integrals are performed with these approximate expressions, and only the 
terms proportional to logarithms of the cutoff parameters are retained. Terms proportional 
to positive powers of the cutoff parameters are set to zero. In order for the method to yield 
reliable results, the cutoff parameters must be kept small, otherwise the approximations 
would not be valid. 

After the two-to-three particle phase space integrals are performed analytically over 
the singular regions and the soft and collinear poles are exposed, the 0(a 2 s ) virtual gluon- 
exchange loop contributions, if any, are added, and all double poles and single poles of soft 
(infrared) origin are verified to cancel, as they should. The remaining collinear singularities 
are factored into parton distribution and fragmentation functions at an appropriate factor- 
ization or fragmentation scale. We work in the MS scheme. One is left with a set of matrix 
elements for effective two-body final-state processes that depend explicitly on ln<5 s and ln<5 c . 
In addition, the non-singular regions of phase space yield a set of three-body final-state 
matrix elements which, when integrated over phase space by Monte Carlo techniques, have 
an implicit dependence on these same logarithms. The signs are such that the dependences 
on In 5 S and In 5 C cancel between the two-body and three-body contributions. The physical 
cross sections are independent of these arbitrary cutoff parameters over wide ranges. In our 
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numerical work, we varied S s and S c over suitable ranges and found quite stable results, as 
is shown in Sec. IV. 

At the level of two-body final-state matrix elements, as in leading-order calculations, it 
is a simple matter to impose selections on kinematic variables similar to those made in ex- 
periments and to calculate different observables. The same is not the case when we consider 
three-body final-state processes. The standard analytic techniques required to obtain dif- 
ferential cross sections of empirical interest often involve complex Jacobian transformations, 
and the phase space integrals can sometimes be done analytically only when specific limits 
of integration are involved. Fully analytic methods of performing calculations for physical 
processes, although in some cases desirable, can be rather restricted in their usefulness when 
it is desirable and sometimes even unavoidable that kinematic selections be made to model 
experimental cuts. The combined analytic and Monte Carlo method is particularly versatile 
in that it provides a means to calculate cross sections differential in many variables at once, 
and to apply cuts on the kinematic variables to match those made in the experiments. The 
phase space integrals are performed numerically after all singularities have been handled 
analytically. 

Our calculation of photon plus charm quark production proceeds along lines similar to 
that for inclusive direct photon production described in Ref. f|, but with a few important 
differences. Since we are interested in observing a final charm quark as well as the photon, 
we cannot integrate analytically over phase space in the limit in which the charm quark is 
collinear to a hard gluon or an anti-charm quark. This situation occurs in the eg, cc, and 
qq initiated processes of Eq. ( |2.1|) . The expression in the appendix of Ref. || for the sum of 
all effective two-body contributions cannot be used in our calculation. We recalculated this 
expression using the three-body matrix elements and virtual gluon-exchange contributions 
from Ref. and we provide the results in our appendix. We also discuss and present in 



the appendix the final-state collinear remnants that give the charm quark momentum distri- 
bution in the limit that the charm quark is produced collinearly with an anti-charm quark 
or a gluon. The singularities in these cases are factored into charm quark fragmentation 
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functions. 



III. CONTRIBUTIONS THROUGH NEXT-TO-LEADING ORDER 

A. Leading order contributions 

In leading order in perturbative QCD, only one direct subprocess contributes to the 
hard-scattering cross section, the QCD Compton process eg —>■ 7c, unlike the case for 
single inclusive prompt photon production, where the annihilation process qq — ► 'jg also 
contributes. Since the leading order direct partonic subprocess has a two-body final state, 
the photon and c quark are produced with balancing transverse momenta. In addition, there 
are effectively leading-order contributions in which the photon is produced by fragmentation 
from a final-state parton. These are 



c + g -> 


# + c 


9 + 9 -> 


c + c 


c + q — > 


c + g 


c + q — > 


c + g 


c + c — > 


c + c 


c + c — > 


c + c 


g + g -> 


c + c. 



(3.1) 



If the photon is to be isolated from the observed charm quark, it arises from fragmentation 
of the gluon g and the non-charm quark q, respectively, in the cases of the first, third and 
fourth processes. In the other cases it is produced by fragmentation of one of the (anti) charm 
quarks. 

In a fully consistent next-to-leading calculation, one should calculate the subprocesses in 
Eq. (|3.1|) to 0(a%), since the photon fragmentation functions that are convoluted with the 
hard subprocess cross sections are of 0(a em /a s ). For simplicity, we include them in 0(a 2 s ) 
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only. In fact, next-to- leading order fragmentation contributions to single prompt photon 



production have been included only once before []13|. We expect the next-to- leading order 
corrections to the fragmentation contributions to be insignificant numerically especially after 
isolation cuts are imposed. 

B. Next-to-leading order contributions 

There are two classes of contributions in next-to-leading order. First there are the virtual 
gluon exchange corrections to the lowest order process, eg — > 7c. Examples are shown in 
Fig. 1(b). These amplitudes interfere with the Born amplitudes and contribute at 0{a em Q 2 s ). 
They were calculated twice before |]8|JT3|1 . At next-to- leading order there are also three-body 
final-state contributions, listed in Eq. Q2.1J) . The matrix elements for these are also taken 



from Ref. [T3|], where they are calculated for single inclusive prompt photon production. 

The main task of our calculation is to integrate the three-body matrix elements over the 
phase space of the unobserved particle in the final state. The situation here is different from 
the standard case of single inclusive particle production because we wish to retain as much 
control as possible over the kinematic variables of a second particle in the final state, while 
at the same time integrating over enough of the phase space to ensure cancellation of all 
infrared and collinear divergences, inherent when massless particles are assumed. All the 
processes of Eq. (|2.1| ), except the first, involve collinear singularities but no soft singularities. 
These collinear singularities must be exposed and factored as explained in Sec. II. The results 
of these calculations are listed in the appendix. 

At 0{af) there are, in addition, fragmentation processes in which the hard-scattering 



two-particle final-state subprocesses 



c + 9 ~ > 7 + c 
c + c — > 7 + g 

q + q^i + g (3.2) 



s 



are followed by fragmentation processes c — > cX, in the case of the first subprocess, and 
g — ► cX in the cases of the last two. These should be included because we have factored 
the collinear singularities in the corresponding three-body final-state processes into non- 
perturbative fragmentation functions for production of a charm quark from a particular 
parton. As a first approximation, we estimate these fragmentation functions by 

D c/C (z,^) = ^p-P qq (z), (3.3) 



and 



D c/g (z,fx 2 ) = ^±P qg (z), (3.4) 



where Pij(z) are the lowest order splitting functions for parton j into parton i |14}| ; and 
a s (fi 2 ) is the strong coupling strength. 



IV. NUMERICAL RESULTS 

In this section we present and discuss several differential cross sections for the joint pro- 
duction of a charm quark and a photon at large values of transverse momentum. All results 
are displayed for pp collisions at the center-of-mass energy = 1.8 TeV appropriate for the 
CDF and DO experiments at Fermilab. To obtain the differential cross sections presented 
in this paper, we convolute our hard-scattering matrix elements with the CTEQ3M parton 
densities ||15|| . Very similar differential distributions may be obtained if other parton sets are 
used instead, with quantitative differences reflecting differences among charm quark densi- 
ties in the different sets |]]. We set the renormalization, factorization, and fragmentation 
scales to a common value fx = pj, in most of our calculations. Dependence on fx is examined 
in one of the figures below. Since there are two particles in the final state, the charm quark 
and the photon, both of whose transverse momenta are large, an alternative choice might 
be fx = Pt or some function of p\ and p^. The results of our calculations show that the 
magnitudes of p? and p c T tend to be comparable and that dependence of the cross sections 
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on \x is slight. Therefore, choices of /x different from // = should not produce significantly 
different answers, and we have verified this supposition in representative cases. 

In addition to showing distributions in the rapidities and transverse momenta of the 
charm quark and the photon, we also discuss distributions in the ratio z, where z is defined 



as 



z = f4 1) 

This ratio z is not to be confused with the variable z in the fragmentation functions, 
Eqs. ( |3.3|) and Q3.4|) . In collider experiments a photon is observed and its momentum is 
well measured only when the photon is isolated from neighboring hadrons. In our calcula- 
tion, we impose isolation in terms of the cone variable R: 



(Ay) 2 + (A0) 2 < R. (4.2) 

In Eq. (|4.2|) , Ay (A0) is the difference between the rapidity (azimuthal angle in the transverse 
plane) of the photon and that of any parton in the final state. The photon is said to be 
isolated in a cone of size R if the ratio of the hadronic energy in the cone and the transverse 
momentum of the photon does not exceed e = 2GeV/pj>- We show distributions for the 
choices R = 0.7 and R = 0.4 typical of current experiments. 

Our first figure in this section, Fig. 2, is an examination of the numerical stability of 
the overall cross section when the cutoff parameters S s and 6 C are varied over appropriate 
ranges. This figure indicates that the combination of analytic and Monte Carlo methods 
yields consistent numerical results for a broad range of the parameters. For the subsequent 
figures, we use S s =0.01, and 6 C =0.001. We return to a brief examination of dependence on 
cutoff parameters when we discuss distributions in the variable z in Fig. 9. 

In Fig. 3 we show the differential cross section as a function of the transverse momentum 
of the charm quark p^, having restricted the transverse momentum of the photon to the 
range 15 < p^ < 45 GeV typical of current hadron collider experiments. The rapidities of 
the charm quark and photon are restricted to the central region, — 1 < y 1,c < 1 in order 
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to mimic the central region coverage of major collider detectors. The solid curve shows 
our prediction when no further selections are made other than those mentioned just above. 
Distributions are presented in the figure for various selections on other kinematic variables. 
For the dashed curve, the variable z of Eq. (|4.1| ) is restricted to z > 0.1. This selection 
on z places the photon and charm quark in opposite hemispheres and results in a modest 
reduction in overall rate. Retaining this cut on z, we examine the effects of isolation of the 
photon and obtain the results shown by the dotted and dot-dashed curves, for cone sizes of 
R — 0.4 and 0.7 respectively. 

A common and notable feature of the curves in Fig. 3 is that slopes change near p T = 15 
and 45 GeV. There is a simple reason for this behavior. The contributions to the cross section 
from two-particle final states produce kinematic configurations in which the photon and 
charm quark have values of px that are equal in magnitude but opposite in sign. Therefore 
two-body final-state processes cannot contribute in the regions p T < 15 GeV and p? > 
45 GeV. Only the three-body final-state processes contribute to the cross section in these 
regions. The steeper fall of the cross section in either direction away from the region 15 < 
Pt < 45 GeV reflects the decreasing likelihood that the photon and charm quark have 
substantially different values of transverse momentum. 

Another feature of the results shown in Fig. 3 is that the effect of isolation diminishes 
as p T is decreased. Isolation affects the cross section principally when the third parton in 
final state enters the photon isolation cone and carries transverse momentum greater than 
the energy resolution threshold. With p T fixed above 15 GeV, the third parton must be in 
the charm quark's hemisphere when p T is small in order to balance p T . When p T is large, 
the third parton is free to enter the photon isolation cone. 

In Fig. 4, we show the cross section differential in p T . The cuts made are the same 
as those for Fig. 3, but in this case the charm quark's transverse momentum is restricted 
between 15 and 45 GeV. The explanation for the change in behavior of the distributions 
above p T = 45 GeV and below p T = 15 GeV is the same as for Fig. 3. An obvious difference 
between Figs. 3 and 4 is that the effect of photon isolation is most significant in the region of 
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small ]P T . The explanation is, again, that with p 1 ^ restricted above 15 GeV, the third parton 
in the final state will be found in the photon hemisphere when pj, is small. It is likely to be 
in the photon isolation cone, and the configuration will be rejected by the isolation cuts. 

To examine further the effects of selections on the charm quark momentum, we present in 
Fig. 5 the cross section differential in the transverse momentum of the photon for a different 
set of cuts. The photon's rapidity is limited to the range —0.5 < y 1 < 0.5, and the ratio z 
is restricted to 0.2 < z < 2.0. These cuts are similar to those of our analytic paper Q. In 
Fig. 5, the solid curve represents the cross section with no isolation cuts imposed, and the 
dashed curve shows the isolated cross section. The dot-dashed curve is the leading order 
prediction, with photon isolation imposed. The behavior seen in Fig. 5 is clearly different 
from that of Fig. 4 in that the cross section does not fall off in the region of small p^, as 
expected, since there in no selection in Fig. 5 on p c T (other than the selection on z). 

In Fig. 6 we show the distribution in the rapidity of the charm quark, y c , for different 
cuts on the photon's transverse momentum. In all cases — 1 < y 1 < 1. The distribution 
in y c is fairly broad, with full-width at half-maximum of about 3.2 units in rapidity. The 
dashed and dot-dashed curves show that the distribution in y c may broaden somewhat as 
Pj, is increased. 

The structure of the QCD hard-scattering matrix element produces positive correlations 
in rapidity at collider energies. To examine correlations more precisely, we study the 
cross section as a function of the difference of the rapidities of the photon and charm quark. 
Results are shown in Fig. 7. In Fig. 7(a), we see that the distribution in Ay is narrower 
than the corresponding distribution in y c shown in Fig. 6. The broader distribution in 
Fig. 6 results from a spread of the approximately Gaussian and relatively narrow dynamical 
distribution of Fig. 7(a) over the range — 1 < y 1 < 1. In Fig. 7(b), we select photons whose 
rapidities are in the forward hemisphere, 1.0 < y 1 < 2.0. We observe that the peak in the 
Ay distribution remains close to Ay = 0, reflecting the predicted |16[ positive dynamical 



correlations, but with the typical value of y c lagging somewhat behind that of the selected 



y 1 - 
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In Fig. 8, we display the differential cross section in y c itself, for two intervals of y 1 in 
the forward rapidity region. These distributions show how the typical rapidity of the charm 
quark follows that of the photon. 

The dependence of the cross section on the variable z, defined in Eq. ( |4.1| ), is indicative 
of the imbalance in transverse momentum of the charm quark and the photon. For two- 
body processes, such as the leading-order Compton subprocess gc — > 7c, the photon and 
charm quark have balancing transverse momenta, and the distribution is a 5-function in z, 
5(1 — z). Contributions away from z = 1 are due to the higher order three-body contributions 
or to fragmentation processes. As discussed in Ref. Q, the photon fragmentation processes 
contribute in the region z > 1 only. Processes in which the charm quark is produced via 
fragmentation contribute in the region < z < 1. We thus expect that the effect of photon 
isolation will be observed in the region z > 1. This expectation is confirmed in the results of 
Fig. 9(a) that show the cross section as a function of z for the non-isolated and isolated cases. 
The solid histogram in Fig. 9(a) agrees quantitatively with the corresponding histogram in 
our analytic paper [|IJ except for differences associated with the different choice of parton 
densities. 

The ^-function behavior at leading order is, of course, moderated by non-perturbative 
effects associated with "intrinsic" transverse momentum of the initial partons as well as by 
next-to-leading order perturbative contributions. In this paper, we are working in the usual 
purely perturbative framework in which the initial partons are assumed to be collinear. For a 
three-parton final state, the region of z near unity is the region in which one of the three final 
partons becomes soft. Sensitivity to soft-gluon effects and the necessity for resummation 
procedures is a common limitation when one considers next-to-leading order contributions 
to an observable that is proportional to a 5-function in leading order. In our calculation 
the soft gluon corrections to the three-body processes are considered as effective two-body 
contributions, as discussed in section II. These contribute to the cross section at 2 = 1, 
meaning that all dependence on the soft cutoff parameter S s is concentrated at z — 1. In 
Fig. 9(b) we show the distribution in z for different values of S s . The cross section is fairly 
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independent of this parameter except when it becomes larger than about 0.02. Similarly, 
in Fig. 9(c), we examine dependence on 8 C , the collinear cutoff parameter. We find fair 
stability over a reasonable range of values of 5 C . The observed 5 C variation is similar to 
the scale dependence seen in Fig. 6 of Ref. |1|]. The two have a related physical origin. In 
the approach of Ref. []IJ, the choice of a scale /i partitions the calculation arbitrarily into 
two-body collinear and three-body final-state contributions. If fi is increased, more of the 
next-to-leading order QCD contributions are placed into the parton distributions (collinear 
kinematics), and less into the exact three-body kinematics. Likewise, here the parameter 
S c partitions the overall QCD matrix element into two-body collinear and three-body final- 
state contributions. It therefore is not a surprise to see the same type of variation with S c 
in Fig. 9(c) as is exhibited by scale variation. 

The infrared sensitivity of the distribution in z, reflected in the S s dependence discussed 
above, and in the scale dependence examined in Ref. means that the z distribution at 
may not be calculated sufficiently reliably at next-to-leading order. Resummation of the 
effects of soft gluon radiation are required, particularly in the region near z = 1. 

The renormalization/factorization/fragmentation scale dependence of our isolated cross 
section is illustrated in Fig. 10 as a function of p c T . The cuts are those of Fig. 3 for cone 
size R = 0.7. All three scales are varied simultaneously, fi = np? T . The dependence on 
scale is negligible in the region 15 < p c T < 45 GeV, whereas outside this region we see some 
dependence. As remarked earlier in our discussion of Fig. 3, both two-body and three-body 
processes contribute to the cross section in the region 15 < p 1 ^ < 45. In this region, the 
next-to-leading order process is complete in the sense that there are both 0(a s ) and 0(a 2 s ) 
contributions. The factorization scales in the hard subprocess cross section compensate and 
cancel those in the structure and fragmentation functions. In the regions p c T < 15 GeV and 
Pt > 45 GeV, only three-body processes contribute. There are no factorization scales in the 
hard subprocess cross sections, except for photon fragmentation scales in the latter region, to 
compensate the scale dependence of the structure functions. The absence of compensating 
terms results in the observed scale dependence in these regions, particularly in the region of 
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small pt- 

A useful measure of the importance of next-to-leading order contributions is the U K- 
factor" , defined as the ratio of the full cross section through next-to-leading order to the full 
leading-order cross section, with fragmentation included. We provide values of K in Fig. 
11 as a function of pj- for two different sets of kinematic selections. Both curves represent 
isolated photon cross sections with R = 0.7. The solid curve is the X-factor appropriate 
to the selections of our Fig. 5. A value of K near 1.5 is in agreement with preliminary 
experimental indications ||. The dashed curve in Fig. 11 is the A"-factor for the cross 
section with the cuts specified in Fig. 4. 



V. SUMMARY AND DISCUSSION 

In this paper we present the results of a calculation of the inclusive production of a 
prompt photon in association with a heavy quark at large values of transverse momentum. 
This analysis is done at next-to- leading order in perturbative QCD. We employ a combi- 
nation of analytic and Monte Carlo integration methods in which infrared and collinear 
singularities of the next-to-leading order matrix elements are handled properly. Our re- 
sults agree quantitatively with those we obtained using purely analytic methods , as they 
should, but the combination of analytic and Monte Carlo methods used in this paper is more 
versatile. We provide differential cross sections in transverse momenta and rapidity, includ- 
ing photon isolation restrictions, that should facilitate contact with experimental results at 
hadron collider energies. We show that the study of two-particle inclusive distributions, 
with specification of the momentum variables of both the final prompt photon and the final 



heavy quark, tests correlations inherent in the QCD matrix elements |16| and should provide 
a means for measuring the charm quark density in the nucleon [IJ. Our results are presented 
in terms of the transverse momentum of the charm quark. In a typical experiment 0, the 
momentum of the quark may be inferred from the momentum of prompt lepton decay prod- 
ucts or the momentum of charm mesons, such as D*'s. Alternatively, our distributions in 
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Pt may be convoluted with charm quark fragmentation functions, deduced from, e.g., e + e 
annihilation data, to provide distributions for the prompt leptons or D*'s. 
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APPENDIX A: ANALYTIC EXPRESSIONS 



In order to make this paper reasonably self-contained, we collect in this appendix all 
the formulae we use in the calculation. If we label the momenta for the generic three-body 
final-state process by 

Pi +P2 -> P3 + Pi + P5, (Al) 

where p 3 denotes the photon and p 4 denotes the observed charm quark, we can define the 
Mandelstam invariants 

i = ( Pl - p 3 y 

U = (pa - Ps) 2 

s = (pi+P2) 2 - (A2) 
We express the two-body final-state cross sections in terms of the scaled variable v, where 

v=l+ t -. (A3) 



1. Two-body contributions 

The effective two-body contribution includes the 0(a 2 ) virtual gluon-exchange loop con- 
tributions and the soft and/or collinear parts of the three-body contributions (this remark 
applies to initial-state collinear contributions only, as explained later). After all soft pole 
singularities are cancelled and all collinear pole singularities are facto, the effective two-body 
contribution is expressed as 



&2,bod y (A + B^'j + c + X) = J dvdxidx 2 



j cg->jc 
a<J coU 

dv 



da 



HO 



+ tf( Xl , M 2 )ff(x 2 , M 2 )^-j—(cg - 7 c) 



(A4) 



plus terms in which the beam and target are interchanged. In this section we use subscript 
c to refer to the charm (or anti-charm) quark. 
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We define 



T, 



2-2v + v 2 



eg 



(A5) 



and v i = 1 — v . In Eq. (0\ 



-(eg -> 7c) 



,2 



1 /iVp 

2 (^ T ^ - 14C F T C9 + 2iV c ln^ + 4(2C F + jV c ) 



x T cg In 2 5 S - ^iV>T C9 In ^ + C F T Cfl (3 + 4 In 6.) In ^+ 
C F T C9 (3 + 41n5 s ) In + ^iV c T C9 (ll + 12 In*,) In ^+ 
-(HJVo - 2N F )T cg In 4 + 4(2C F - N c ) Inv + AT cg \n6 s (N c lnv x + 
2C F In u - iV c In u) + N c In 2 v x (l + v) + (2C F - N c ) In 2 v 
(2-2v + 3v 2 ) , o ^ (1 + 2 V ) _ iVc7r 2(z 1 - 2v + 4v2 ) 



x 



+ 2C F lnv 1 - 



+ 2C^ (1 - 4 " + 5 " 2) + 2C fl n 2 , 1 ( " 2 + " ;) 



3f i 



3t>i 

-2(2C F - Ar c )hn;i 



x ln „i^±^ll + 2(2C F - N c )T cg U 2 (l -v) + 2N c T cg U 2 (v) 



(A6) 



The scales M and M" are the factorization and fragmentation scales, respectively, on the 
initial parton and final-state charm quark legs, and fi is the renormalization scale. C F = 4/3 
is the quark-gluon vertex color factor, Nc = 3 is the number of colors, e c is the fractional 
charge of the charm quark, S s and 5 C are the soft and collinear cutoff parameters defined in 
Sec. II, Li 2 (x) is the dilogarithm function, and a em is the electromagnetic coupling constant. 
The remnants of the factorization of the hard collinear singularities are 



d<7coU JC Ole.m.oiiel 1 



dv 



2.S 



X 



T 

N C C9 

r2\ 



M 2 ) ( / ^ ^/f ( ^, M 2 ) P gq (z) + fl dZ < B fX2 1 11 



-f 

J a 



+ /f(*2,M 2 



x 2 Z 



X'! 



Si 



(A7) 



The last distribution, f£(xi, M 2 ), in Eq. ( |A7| ) implies a sum over the flavors of quarks 
from the cc, cq, and cq initial states. The remaining two processes, cc and qq, do not 
have initial state collinear singularities and thus do not contribute to this part of the cross 



section. Equation ( |A7|) differs slightly from that given in Ref. ||, where the upper limit 
of integration was always taken to be 1 — 5 a . In principle, this is incorrect since it implies 
that processes other than the eg initiated process may have soft singularities, but the error 
introduced by that approximation is very small for small 5 S . We use the correct expression 
in our calculation. 

The splitting functions Py, listed in the appendix of Ref. f|, are included here for 
completeness. 



P^z) = P l3 {z) In (^^J - Pljiz). 



(A8) 



The functions Pij{z) are the usual Altarelli-Parisi splitting functions in 4 — 2e dimensions 
and are 

'1 + z 2 



P qq (z, e) = C F 



2(1 -e 



P gg (z,e) =2N C 



<1 - z) 

1 — Z 

z 2 + (1 - z) 
1 - z 



1 — z z 
1 + (1 - zf 

z 



+ z(l-z) 



€Z 



(A9) 



The functions PL(z) are defined by the relation 



P ij (z,e) = P ii (z)+ePlAz) 



(A10) 



2. Pseudo-Two-Body Contributions 

Since we are interested in distributions in the kinematic variables of two final-state 
partons, the photon and the c quark, we can define variables that depend on the momenta 
of both. An example is the variable z, defined in Eq. ( |4.1| ). Whenever there is a third parton 
in the final state, the distribution in z (or in other analogous variables) will differ from a 
delta-function when the third parton carries a finite momentum, even if it is collinear to one 
of the other final partons. For this reason we designate as "pseudo-two-body contributions" 
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those for which the third parton is collinear to either the final photon or the charm quark. 
These contributions are expressed, respectively, by the equations 



"i/cou = £/ f a A (x u M*)f b B (x 2 ,M*) 

abcq 

x —(ab — > cq\dx\dx 2 dzdv ^ 
dv 



P iq (z) In 



2(1 — ,z)5 c 



M' 2 



(All) 



and 



c I coll 



£ I f^(x 1 ,M 2 )f b B (x 2 ,M 2 )P cd (z,M" 2 )^-(ab - 1 d)dx 1 dx 2 dzdv. (A12) 



The functions P 1 / q {z) and P^i q (z) are the quark-to-photon splitting function and 0(e) piece, 
respectively. They have the same form as P gq , with the color factor replaced by the square 
of the quark charge. The scale M' is the fragmentation scale for quark fragmentation into 
a photon. 



P cd (z } M" 2 ) =P cd (z) In 



z(l — z)5 c 



M" 2 



(A13) 



where P c d(z) represents the splitting functions P qg (z) and P qq (z) of Eq. ( |A9| ) along with the 
'primed' pieces. In Eq. ( |A11| ), q can be a charm or anti-charm quark, or a (anti) quark of 
any flavor in the case of cq — > cq. In Eq. ( |A12| ), d is either a gluon or (anti) charm quark. 

These contributions are usually referred to as the remnants of the factorization of the 
hard collinear singularities and are regarded as two-body processes, or as parts of the frag- 
mentation contributions because of their dependence on the factorization scales. We prefer 
to regard them as pseudo-two-body contributions. When we examine either the charm quark 
or photon momentum distributions, these contributions populate the same regions of phase 
space as the other three-body contributions in Eq. (|A17|) , unlike the effective two-body con- 
tributions. The pseudo-two-body contributions are usually negative in overall sign due to 
the large logarithms of the cut-off parameters, as are the two-body contributions discussed 
above. 
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3. Photon Fragmentation Contributions 



As mentioned in Sec. Ill, we include the quark-to-photon and gluon-to-photon frag- 
mentation contributions at leading order only. We convolute the 2^2 hard scattering 
subprocess cross sections for the processes listed in Eq. ( |3.1| ) with photon fragmentation 
functions D 1 u(z, M' 2 ); M' is the fragmentation scale, the same scale at which we subtract 
the collinear singularities on the photon leg of the three-body processes. The expression for 
the cross section is 

°fra 9 h = E / fa(*U M 2 )f?{x 2 , M 2 )D^(z, M' 2 )^{ab -> jc)dx 1 dx 2 dzdv . (A14) 

abj J dV 

The matrix elements for the hard subprocess cross section are well known and can be found, 
for example, in Ref. f4[. 



4. Charm Fragmentation Contributions 

In integrating some of the three-body matrix elements over phase space we encounter 
configurations in which the charm quark is produced collinearly with an anti-charm quark 
or a gluon in the final state. These situations lead to a collinear singularity in the massless 
approximation. They occur for the processes of Eq. ( [2. lap , ( 2.1c ), and (|2.H]) . We factor 
these singularities into a fragmentation function D c /i(z, M" 2 ) for parton i to produce a charm 
quark with momentum fraction z. The contributing subprocess cross sections are 



• . 2C F nctat s e 2 ( v 1 — v\ 

-{qq -> 7# = — H - - + 

Nn s \ 1 — V v J 



da ^ N 2C F Ttaa s e z q ( v \ 

dv i\c 
da. . iraa s e 2 (l + (1 - v) 2 \ 

-r(qg - 7?) = — 7 • ( A15 ) 



dv Ncs \ 1 — v 

The physical cross section is given by 

da 



°fra 9 /c = E / fa(*u M 2 )f^(x 2 , M 2 )D c/d (z, M" 2 )^(ab - jd)dx l dx 2 dzdv . (A16) 

abd 
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5. Three-body contributions 



The non-collinear three-body final-state contributions are calculated from the expression 
o-i-body = E / fa(xu M 2 )f b B (x 2 , M 2 )da(ab -> rycd)dx 1 dx 2 dn, (A17) 



abd 



with Q representing the angles and other variables that are integrated over. Whenever an 
invariant sy or t{j falls into a collinear or soft region of phase space, that contribution from 
the subprocess is excluded. The three-body contribution shows no dependence on the fac- 
torization scale of the final-state charm or photon legs, although we have factored collinear 
singularities at scales M" and M', respectively, on these legs of the three-body subprocesses. 
However, Eq. (|A17 ) does contain implicit logarithmic dependence on the soft and, in partic- 



ular, the collinear cutoffs discussed in Sec. II. Both collinear cutoff and factorization scale 
dependences are contained in the pseudo-two-body contributions discussed above. 
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FIGURE CAPTIONS 

[1] (a) Lowest order Feynman diagrams for 7 plus c quark production; k\ and ki 
are the four-vector momenta of the photon and charm quark, (b) Examples of 
virtual corrections to the lowest order diagrams, (c) Examples of next-to-leading 
order three-body final-state diagrams for the gc initial state. 

[2] Study of the dependence of the cross section on the Monte Carlo cutoff param- 
eters 5 S and 5 C . Shown is the cross section for p + p — > 7 + c + X at y/s = 1.8 
TeV with the transverse momenta of the photon and charm quark restricted to 
the the interval 10 < p T < 50 GeV and the rapidities of the photon and charm 
quark restricted to the interval —3.0 < y < 3.0. 

[3] Cross section da/dp^ as a function of the transverse momentum of the charm 
quark for p + p — » 7 + c + X at y/s = 1.8 TeV. The transverse momentum of 
the photon is restricted to the interval 15 < p? < 45 GeV, and the rapidities 
of the photon and charm quark are restricted to the interval —1.0 < y < 1.0. 
Four curves are drawn. The solid curve shows the cross section with no further 
restrictions. The dashed curve indicates the result after the additional selection 
is made that z > 0.1; the ratio z is defined in the text. The dotted and dash- 
dot curves display the results after photon isolation restrictions are applied, in 
addition to the cut on z. 

[4] Cross section da/dp? as a function of the transverse momentum of the photon 
for p + p — > 7 + c + X at =1-8 TeV. The transverse momentum of the 
charm quark is restricted to the interval 15 < p c T < 45 GeV, and the rapidities 
of the photon and charm quark are restricted to the interval —1.0 < y < 1.0. 
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Four curves are drawn. The solid curve shows the cross section with no further 
restrictions. The dashed curve indicates the result after the additional selection 
is made that z > 0.1; the ratio z is defined in the text. The dotted and dash- 
dot curves display the results after photon isolation restrictions are applied, in 
addition to the cut on z. 

[5] Cross section da /dp], as a function of the transverse momentum of the photon 
for p + p — >7 + c + X at ^/s = 1.8 TeV. The transverse momentum and rapidity 
of the charm quark are not restricted, but the rapidity of the photon is limited 
to the interval —0.5 < y < 0.5. The ratio z is restricted to the interval 0.2 < z < 
2.0. Three curves are drawn. The solid curve shows the cross section with no 
further restrictions. The dashed curve indicates the result after photon isolation 
is imposed, with R = 0.7, and the dot-dashed curve is the leading order cross 
section with photon isolation imposed, R = 0.7. 

[6] Cross section da/ dy c as a function of the rapidity of the charm quark for p + p — > 
7 + c + X at y/s =1.8 TeV. The photon rapidity is restricted to —1.0 < y 1 < 1.0, 
and photon isolation is imposed, with R = 0.7. There is no restriction on p^, 
but the ratio z is restricted to z > 0.1. Curves are shown for three selections 
on the transverse momentum of the photon. The solid, dashed, and dot-dashed 
curves correspond to the selections 15 < p\ < 45 GeV, 15 < < 25GeV, 
and 35 < p\ < 45 GeV, respectively. For ease of comparison of shapes, the 
dot-dashed curve has been multiplied by 10. 

[7] Cross section da/dAy as a function of the difference Ay = y 1 —y c of the rapidities 
of the photon and charm quark, for p + p — > 7 + c + X at =1.8 TeV. The 
ratio z is restricted to z > 0.1, and photon isolation is imposed, with R = 
0.7. The transverse momentum of the photon is selected to be in the interval 
15 < < 45 GeV. In (a), the photon rapidity is restricted to —1.0 < y 1 < 1.0; 
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in (b) 1.0 < y 1 < 2.0. In (b), the dashed curve shows the behavior at leading 
order. 

[8] Cross section da/ dy c as a function of the rapidity of the charm quark for p + p — > 
7 + c + X at y/s = 1.8 TeV. Photon isolation is imposed, with R = 0.7; z > 0.1; 
and 15 < p? < 45 GeV. The solid curve shows the result when the photon 
rapidity is restricted to 1.0 < y"' < 2.0, and the dashed curve displays the result 
for 2.0 < y' < 3.0. 

[9] Cross section da jdp^dy'dz as a function of the ratio zioip + p^^ + c + X 
at y/s = 1.8 TeV. The transverse momentum and rapidity of the photon are 
averaged over the intervals 14 < p\ < 16 GeV and —0.5 < y 1 < 0.5. In (a), we 
illustrate the effects of photon isolation by comparing the distributions in z with 
and without the isolation restriction. The solid histograms in (b) and (c) show 
the results of our calculation for our standard Monte Carlo integration cutoff 
parameters 5 S = 0.01 and 5 C = 0.001. In (b) and (c), photon isolation is imposed 
with R = 0.7. In (b), we display the dependence of the final cross section on the 
Monte Carlo integration cutoff parameter S s , having fixed 5 C = 0.001. In (c), we 
show the dependence of the final cross section on the cutoff parameter S c , having 
fixed 5 S = 0.01. 

[10] The renormalization/factorization scale (/i) dependence of the cross section is 
displayed. Shown is the transverse momentum dependence of da /dp^ for three 
values of fx/p?'- 0.5, 1.0, and 2. The transverse momentum of the photon is 
restricted to the interval 15 < p? < 45 GeV, and the rapidities of the photon 
and charm quark are restricted to the interval —1.0 < y < 1.0. 

[11] The K factor defined in the text is shown as a function of pj.. Photons are 
isolated with R = 0.7. The solid line corresponds to selections analogous to 
those of our analytic paper [Ref. 1]: —0.5 < y 1 < 0.5 and 0.2 < z < 2.0, 
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but with no restrictions on y c or p?. The dashed line shows how results change 
when the transverse momentum of the charm quark is restricted to the interval 
15 < < 45 GeV with, in addition, z > 0.1 and the rapidities of the photon 
and charm quark limited to the interval —1.0 < y < 1.0. 
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